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Abstract 

The abundances of the chemical elements observed at the surface 
of metal-poor stars are not always representative of their initial values. 
During stellar evolution, various physical processes modify their internal 
composition. In this short paper, in honor of George W. Preston, I 
remind the importance of atomic diffusion, and the possible effects of 
accretion processes which may lead to double-diffusive (thermohaline) 
instabilities. I discuss the consequences of these processes and compare 
them to current observations. 

1. Introduction 

The chemical abundances observed in metal poor stars are often considered as a 
clue for the understanding of the chemical evolution of galaxies, and they are generally 
assumed to be representative of their initial values. This assumption does not take 
into account the basics of stellar physics, which lead to the conclusion that the element 
abundances must change with time during stellar evolution, due to atomic diffusion 
acting alone or combined with the various macroscopic processes which may occur 
inside the stars. 

An interesting case is that of lithium in halo stars. Stellar physicists claimed for a 
long time that the abundances of lithium in these stars cannot be the original ones, 
the computations showing that it has decreased by at least a factor two sin ce th e 
origin, du e to a tomic diffusion coupled wit h nuc lear re action s (Michaud et al. | (jl984 ). 



Vauclair | fll988f ). Vauclair & Charbonnel| f1l995f ) and| fll998f )l In spite of these results 



the observed lithium value was still presented by the observers as the primordial one, 
and used as a constraint for the Big Bang nucleosynthesis. Later on, the space mission 
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WMAP proved that the baryonic density obtained from the cosmological background is 
different from that obtained with the observed lithium value, and that it is compatible 
with a smaller primordial value, as predicted by stellar physicists. 

Atomic diffusion is more complex to compute for heavy elements than for light 
ones, because of their more complicated atomic structure. It should, however, always 
been taken into account in the interpretation of the observations. In the following, I 
first remind the basics of atomic diffusion, and I give examples of its importance for 
metal-poor stars (Section 2). In section 3, I discuss another process acting on some 
metal poor stars: accretion followed by thermohaline convection. These processes are 
connected, as the abundance gradients induced by atomic diffusion play an important 
role in the extent of the mixing. I give my conclusions in Section 4. Without George 
Preston, I would not have been involved in the problem of accretion onto metal poor 
stars, which is a very interesting subject from both observational and physical view 
points. Much thanks to you, George, on your birthday. 



2. Atomic diffusion in stars 

The fundamental importance of atomic diffus ion inside sta rs was recognized by the 
pioneers of stellar physics like Eddington| ( 1916 ). Chapman | ( 1917 ) already discussed 



the "relative importance of convection and diffusion within a star". At that time, he 
only introduced in the computations the effects of gravitational and thermal diffusion. 
In this case, he found that no heavy elements should be left at the surface of stars, 
which was obviously wrong. 

A first approach of the possible effects of selective radiation pressure, pushing the 
e lemen ts upwards in a way related to their atomic structure, may be found in Eddington 
| ( ]l926 ). He explained how the radiation force is much stronger on heavy atoms than 



on light ones because of line absorption: the heavy elements should be pushed up at 
the surface of stars, with an effect increasing with increasing stellar mass. But this 
was not observed at that time and was put at the back of the mind of astronomers for 
a long time. 

Aller & Chapman | fll960h came back to this problem and computed atomic diffusion 



time scales for metals at the bottom of the solar convective zone. They found that this 
time scale was longer than the solar age, so that there should be no important effect 
of atomic diffusion inside the Sun. 

Some years later, the presence of chemical peculiarities in main sequence stars 
became more and more evident, and the explanation in ter ms of atomic diffusion more 
and more successful (see references in Vauclair & Vauclair J (jl982h ). It was still believed 
however that atomic diffusion occurred only in specific stars, and was suppressed by 
convection or other macroscopic motions in most of the other stars. In particular, 
the small effect of helium and metal diffusion in the Sun was considered as negligible, 
until helioseismology proved in a spectacular way that it really occurs as predicted and 
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has a measurab le infl uence on the internal solar sound velocity (Gough et al. | ( 1996 ). 
Richard et al. I (Il996j )). 

From then on, atomic diffusion in stars has been much studied. The most difficult 
problem in the computations is the evaluation of selective radiative accelerations on 
each ion of each element. This is related to the computations of opacities, as radia- 
tive accelerations act through the absorption of the radiation flux by the ions. At 
the present time, three stellar evolution codes are able to compute ato mic dif fusion 
including radiative accelerations: the Montreal code (C.f. Richard et al. I (20021 )). the 



TGEC code (Toulouse Geneva Evolut ion C ode, C.f. Theado et al. I ( ]201ol ) ) and the 
Yale code (Delahaye & Pinsonneault | (j2005l )). 

Atomic diffusion cannot occur inside convective zones, due to the efficient mixing, 
and it is subject to competition with extra-mixing processes inside radiative zones. 
These processes can be due to rotation-induced mixing, internal waves, mass loss, 
thermohaline convection, etc. Th ey ha ve been extensiv ely d iscussed in the recent 
literature (see Thea do & Vauclair | (l2003h . Thea do et al. I feoiflh . Talon & Charbonnel 



I (l2005h and | (l2008f ). Dennissenkov et alj fcoioh . etc.). Due to all these macroscopic 



processes, atomic diffusion is slowed down and the final results modified, but it is 
generally not completely suppressed. Atomic diffusion alone is parameter free, which 
is not the case for mixing. The detailed observation of the abundances of many elements 
can lead to strong constraints on these macroscopic processes. 

Atomic diffusion in metal poor stars is important for various reasons. 

• The helium settling leads to modifications of the stellar evolution trac ks, an d this 
may change the evaluated age of globular clusters. Van den Berg et alj (12002 ) have 



found that introducing atomic diffusion can reduce the age by about ten percent. 

The presence of diffusion-induced helium gradients below the convective zones can 
have important secondary effects, for example in the case of metal-rich accretion, 
as discussed below (Section 3). 

Beyond the lithium problem, discussed in the introduction, atomic diffusion mod- 
i fies th e det ailed a bundances of heavy elem ents in halo stars (e.g. Richard et al. 
d2002h and I (l2005f ). Korn et al. | (120061 ) and| fj2009h ). Concerning lithium, diffusion 



alone would lead to results incompatible with the observations, as it should be 
more deficient in more massive stars, which is not observed. Mixing processes are 
needed to slow down the selective process. However, such a mixing connects the 
bottom of the convective zone to the lithium nuclear destruction layers, so that in 
any case lithium is depleted. The fact that lithium is not completely absent from 
the spectra of metal deficient stars leads to constraints on the mixing processes 
which must not be too strong. As a consequence, atomic diffusion still plays a 
role and this is r eflect ed on the abundan ces of metals (e.g. discussions in Vauclair 
& Charbonnel I (Il998h and Korn et al. I (120061 )) 
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3. Accretion and thermohaline convection 

The abundances at the surface of stars may sometimes be altered by the accretion 
of external matter. Two interesting cases have recently been studied. The firs t case 
concerns the accretion of planetar y mat ter onto exoplanet-host sta rs. Va uclair I ( 120041 ) 



a nd m ore recently Theado et al. I (l2010h and Theado & Vauclair and Garaud 



I ( 201lh have discussed how the accreted matter is mixed down due to double-diffusive 



(thermohaline) convection induced by the inverse /i-gradient. 

Here I will discuss another case, in which George Preston is particularly interested: 
that of Carbon Enriched Metal Poor Stars (CEMP). 

CEMP stars are metal poor stars ([Fe/H] smaller than -2), with a large overabun- 
dance of carbon with respect to the other heavy elements ([C/Fe] larger tha n +1.0 ). 
The number of these stars increase with decreasing metallicity (Lucatello et al. | (j2006h ) . 
Most of them are rich in s- process elements, forming the so-called CEMP-s group. The 
commonly accepted explanation of the C and s process elements enrichment is these 
stars is that they were submitted to a shower of matter coming from an AGB compan- 
ion which later on can have become a white dwarf. 

A complication comes from the fact that about half of the CEMP-s stars also 
show a high r-process enrichment, incompatible with a pure s-process contribution. 
The nucleosynthesis of the r-process does not come from AGBs, but it is believed to 
come from massive stars exploding as Type II Supernovae. An interpretation of these 
observations is that the r-process enrichment was already presen t in the molecular 
cloud from which the bina ry system formed (Bisterzo et al. I ( 2010l )). 



Stancliffe et al. | ( 2007 ) rejected the accretion hypothesis on the basis that in the 



case of accretion of metal-rich matter, thermohaline convection develops so th at the 
accreted matter is rapidly diluted throughout the star. Thompson et al. I (Hoosl) 



pointed out that the main sequence star which suffers accretion was subject to helium 
settling since its birth so that it could have built a stable /i-gradient in its outer 
layers. Such a gravitational-settling induced /i-gradient can stop at least some of the 
accretion-induced thermohaline mixing. In this case the chemical composition of the 
stellar atmosphere can be modified by the accreted material. Later on, Stancliffe et al. 



(120081 ) performed new computations including gravitational settling and confirmed that 



it may slow down or even inhibit the thermohaline mixin g. In their computations, they 
used the prescription suggested by Charbonnel & Zahn| (12007 ) for Red Giants, which is 



now believed to be overestimated. We can thus infer that the mixing effect in CEMPs 
is still smaller than previously thought, so that the accretion- induced modification of 
the abundances may be larger than previously thought. Let me discuss this in more 
detail 

Thermohaline convection is a well known process in oceanography : warm salted 
layers on the top of cool unsalted ones rapidly diffuse downwards even in the presence 
of stabilizing temperature gradients, due to the different diffusivities of heat and salt. 
When a warm salted blob falls down in cool fresh water, the heat diffuses out more 
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quickly than the salt. The blob goes on falling due to its weight until it mixes with 
the surroundings. This leads to the so-called "salt fingers" . Thermohaline convection 
is a "double diffusive convection". When the gradients are reversed, another kind of 
double diffusive convection occurs, which is generally referred to as "semi convection" . 
The condition for the salt fingers to develop is related to the density variations induced 
by temperature and salinity perturbations. 

Thermohaline convection occurs in stellar radiative zones in the presence of inverse 
/i-gradients. The medium can become dynamically unstable if (Ledoux criterion): 



crit 



V 



ad 



V < 



(1) 



where (j) — {d In p/d In p,) and 6 — (d In p/d In T). When V C ru vanishes, marginal 
stability is achieved and thermohaline convection may begin as a "secular process", 
namely on a thermal time scale (short compared to the stellar lifetime). 

The effects of thermohaline convection as a mixing proce ss in stars are far from 
trivial. First p arametrization recipes were given by Ulrich I (Il972h and Kippenhahn 
et al. | (jl98Cll ). These authors treated thermohaline mixing as a diffusion process 
with a diffusion coefficient depending on the relative size of the "fingers", and more 
specifically on the ratio of their horizontal to vertical scale, which was totally unknown. 
Furt hermo re, a multiplying coefficient was intro duced, of order 12 for Ki ppenh ahn et 
al. | (Il980h . and much larger for Ulrich | (1l972f ). Charbonnel & Zahn | (j2007h found 
that, to account for the abundances observed in red giants, this coefficient should be 
of order 1000. 

More recently, numerical simulations of thermohaline convection i n con ditions closer 
t o the stellar ones than before have been performed by Traxler et al. | (12011 ) and Garaud 
They derive a new diffusion co efficie nt which is of the same order of magnitude 
as found from the Kippen hahn et al. | (198(3) prescription, much below that suggested 
by Charbonnel & Zahn| (l2007h . This result is confi rmed through different numerical 
computati ons b y Denissenkov & Merryfield | (j2011 ). It is also shown by Theado & 
Vauclair | (12011 ) that the large mixing coefficient needed for the Red Giants would lead 
to more lithium depletion than observed in case of planetary material accretion onto 
the star. 



4. Conclusion 

I first met George W. Preston in 1978, when, during a long stay in Cal Tech, I had 
the privilege to do some observations with my husband Gerard Vauclair at the 200 
inch telescope in Mount Wilson. We did spectroscopic observations of Am stars, and 
George was our (extremely nice and efficient) support astronomer. We had afterward 
some festive evenings, doing music together (I still have the american folksong book 
he gave me at that time and I use it for my grand daughters!). At that time, George 
already told us about metal-poor stars, but we were more involved with metal-rich 
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ones... 

Time has passed, and it was a great pleasure for me to be sought by him a few years 
ago about the efficiency of thermohaline convection in CEMPs stars, in case of AGB 
accretion. I remember he really wanted me to discuss about mixing coefficients and time 
scales. I thought about it and then came in my mind the idea of the helium-induced 
stabilizing /z-gradient. Such a situation does not occur for the accretion of planetary 
matter onto exoplanet-host stars, which I studied before, because the stars are very 
young when they suffer the metal-rich accretion. But when old main-sequence stars 
accrete AGB matter, they obviously had time to built such a gradient. Helioseismology 
has proved that helium settling is at work, even in stars with deep convective zone as 
in the Sun. 

This was for me a very interesting subject, which I would not have entered without 
George. The story is not finished, and much work still has to be done on the simulations 
of various situations were accretion, thermohaline mixing and gravitational settling 
work together, besides nuclear reactions, to lead to the rich variety of observed element 
abundances. Nature is complex and fascinating, and we are lucky to work on a better 
understanding of all this complexity. Thanks and long live! 
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